On solvability and integrability of the Rabi model 
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Quasi-exactly solvable Rabi model is investigated within the framework of the Bargmann Hilbert 
space of analytic functions B. On applying the theory of orthogonal polynomials, the eigenvalue 
equation and eigenfunctions are shown to be determined in terms of three systems of monic ortho- 
gonal polynomials. The formal Schweber quantization criterion for an energy variable x, originally 
expressed in terms of infinite continued fractions, can be recast in terms of a meromorphic function 
F{z) = flo + X^J^i Mk/{z — ^k) in the complex plane C with real simple poles and positive residues 
Mk- The zeros of F{x) on the real axis determine the spectrum of the Rabi model. One obtains 
at once that, on the real axis, (i) F(x) monotonically decreases from +00 to —00 between any two 
of its subsequent poles ^k and ^fc+i, (ii) there is exactly one zero of F{x) for x £ i^k,^k+i), and 
(iii) the spectrum corresponding to the zeros of F{x) is necessarily nondegenerate in each of the two 
parity eigenspaces B±. Thereby the calculation of spectra is greatly facilitated. Our results allow 
us to critically examine recent claims regarding solvability and integrability of the Rabi model. 



I. INTRODUCTION 



The Rabi model [2] describes the simplest interaction 
between a cavity mode with a bare frequency w and a 
two-level system with a bare resonance frequency wq. 
The model is characterized by the Hamiltonian 



H = hwa'a + hgai{a) + a) + //era, 



(1) 



where a and are the conventional boson annihilation 
and creation operators satisfying commutation relation 
[a, a^] — 1, g \s a coupling constant, and /i — /iwo/2. 
In what follows we assume the standard representation 
of the Pauli matrices (jj and set the reduced Planck 
constant ft = 1. For dimensionless coupling strength 
K = g/uj < 10~^, the physics of the Rabi model is 
well captured by the analytically solvable approximate 
Jaynes and Cummings (JC) model H, The latter 
is obtained from the former upon applying the rotating 
wave approximation (RWA), whereby the coupling term 
(Ji (a^ + a) in Eq. ^ is replaced by {a+a + a-o,^), where 
a± = ((Ti ± ia-i)!^. Nowdays, solid-state semiconduc- 
tor and superconductor systems have allowed 
the advent of the ultrastrong coupling regime, where the 
dimensionless coupling strength k > 0.1 (IQ] . In this 
regime, the validity of the RWA breaks down and the 
relevant physics can only be described by the full Rabi 
model J|. With new experiments rapidly approaching 
the limit of the deep strong coupling regime characterized 
in that k > 1 pd| . i.e., an order of magnitude stronger 
coupling, the relevance of the Rabi model [l| becomes 
even more prominent. There is every reason to believe 
that ultrastrong and deep strong coupling systems could 
open up a rich vein of research on truly quantum effects 
with implications for quantum information science and 
fundamental quantum optics [6] . 

The Rabi model applies to a great variety of physical 
systems, including cavity and circuit quantum electro- 
dynamics, quantum dots, polaronic physics and trapped 
ions. In spite of recent claims 0, [ill, the model is not 
exactly solvable. Rather it is a typical example of quasi- 




exactly solvable (QES) models in quantum mechanics 
[l3l - [l7t . The QES models are distinguished by the fact 
that a finite number of their eigenvalues and correspond- 
ing eigenfunctions can be determined algebraically (Tsl - 
Il6l | . That is also the case of the Rabi model ^3] ■ Certain 
energy levels of the Rabi model, known as Juddian ex- 
act isolated solutions [l^, can be analytically computed 
42(3l| . whereas the remaining part of the spectrum not 
|20|. Depending on model parameters, the spectrum 
can only be approximated (sometime rather accurately 
- cf. Eqs. (18), (20) and Fig. 3 of Ref. [211]; Eq. (20) 
and Figs. 1,2 of Ref. [l^]). Therefore, any kind of exact 
results involving the Rabi model continues to be of great 
theoretical and experimental value. 

In our earlier work [23l | we studied the Rabi model as 
a member of a more general class TZ of quantum mod- 
els. Models of the class TZ were characterized in that the 
eigenvalue equation 



(2) 



where H denotes a corresponding Hamiltonian, reduces 
in the Bargmann Hilbert space of analytic functions B 
0, [2^ to a three-term difference equation 



■ an(t>n 







(n > 0). (3) 



Here {0n}^o ^^^^ ^^'^ sought expansion coefficients of a 
physical state described by an entire function 



(4) 



Models of the class TZ were then characterized in that ad- 
ditionally the recurrence coefficients have an asymptotic 
powerlike dependence [1^ 



an 



bn" {n — ^ 00), 



(5) 



where a and b are proportionality constants and the ex- 
ponents satisfy 2i; > u and t = <; — v > 1/2 f23l|. The 



2 



spectrum of any quantum model of TZ can be obtained 
as zeros of the transcendental function of a dimensionless 
energy parameter x [IE [1^ , 



F{x) 



■ • Pk, 



(6) 



fe=i 



where F{x) is defined solely in terms of the coefficients 
of the three-term recurrence uM, : 



Pi 



ui = 



bi 

, Pi 

ai 



ui 



1 



1 - ui-ihi / {am^iY 



1, Ul = 1, 



I > 2. 



(7) 



The function F{x) yields both analytic and efficient nu- 
merical representation of the formal Schweber quanti- 
zation criterion expressed in terms of infinite continued 
fractions (cf. Eq. (A.16) of Ref. [2]), 




F{x) 



ai- 02- as- 



(8) 



The Rabi Hamiltonian Hji is known to possess a dis- 
crete Z2-symmetry correspondin g to the constant of mo- 
tion, or parity, fl = exp(i7rJ) where 



J = a^a + i (1 + erg) 



(9) 



is the familiar operator known to generate a continuous 
U{1) symmetry of the JC model [MQ- The Rabi model 
can be characterized by a pair of the three-term recur- 
rences (Eq. (37) of Ref. ^) 



1 



K{n + 1) 



6±(-l)"A]</.± 

1 



1 



0. 



(10) 



Here e = /w, k = g/uj reflects the coupling strength, 
A = n/uj, and ± applies to the positive and negative 
parity states, respectively, with power series expansions 
$±(2) = yi ^n j iiz^ for the respective parity states be- 
ing assumed [23|. The Bargmann Hilbert space of an- 
alytic functions can be thus written as a direct sum 
B ^ ® B- of the parity eigenspaces. The case of a 
displaced harmonic oscillator, which is the exactly solv- 
able limit of for /i = 0, results in A = 0, whereby the 
recurrence PU)) reduces to Eq. (A. 17) of Ref. [2]. 

In the case of the Rabi model, and its special case of 
the displaced harmonic oscillator, it was observed that 
the plots of F{e) corresponding to ([TU)) displayed a se- 
ries of discontinuous branches monotonically decreasing 
between +00 and —00 (see Fig. [T]and Figs. 1,2 of Ref. 
[23j). In the present work the latter property will be 
proven. First we show that F{x), considered as a func- 
tion oi X = e/ n = /g, can be alternatively expressed 



Figure 1. F{x) corresponding to the recurrence (|10p in the 
deep strong coupling regime for k — 1.4, A = 0.4, lo = 1, and 
X = e/K = E/g. 



as the limit of rational functions 



F{x) = lim Fn{x) = ao + lim 



(11) 



where {P„(x)}^=o and {Pi'^(x)}^= Q are associated sys- 
tems of monic orthogonal polynomials. (Monic means 
here that the coefficient of the highest power of x is one.) 
The polynomials of each orthogonal polynomial system 
(OPS) {Pn{x)}^=i and {Pk^\x)}^^, 

• have real and simple zeros, and 

• the zeros of P,jL\(a;) and Pn{x) are interlaced. 

Specifically, denote the zeros of Pn{x) with degree P„ = n 
by Xni < Xn2 < . ■ . < Xnn and the zeros of pj^^}i{x) with 



degree n - 1 by x^^^l^ ^ < x^'^^ 2 < ■ • ■ < <':^i,„_i- Then 
for any fc = 1, 2, . . . , n — 1 



.(1) 



.(1) 



(a) (a) (a) 

■^nk "^n—l.k *^ •^n.k-\-li 



(12) 
(13) 



where a = 0, 1 (for the sake of notation the superscript 
(0) for a = will be suppressed in what follows). For 
each fixed fc, {a;„fc},^j. is a decreasing sequence and the 
limit 

= lim Xnk (14) 

n— >oo 

exists. Additionally, for any finite n the ratio in (|lll) . 
also known as a convergent, enables a partial fraction 
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decomposition (PFD) 



E 

fc=i 



(15) 



The numbers M„fe are all positive^ Mnk > 0, and satisfy 
the condition 



(16) 



fc=i 



Each number Mnk can be shown to correspond to the 
weight corresponding to the zero Xnk in the Gauss 
quadrature formula for the positive definite moment 
functional C associated with the OPS {Pn{x)}^^Q. In 
the case of the displaced harmonic oscillator and the Rabi 
model, 



(n+1) 



Pn+l{Xnk)Pn{Xnk) 

-1 




> 0. 



(17) 



From (|15l) one finds immediately that whenever the 
derivative dFn{x)/dx exists, then 



dFnjx) 
dx 



< 0. 



(18) 



Consequently, between any two subsequent Xnk < Xn,k+i, 
where Fn{x) decreases from +oo to — oo, there is exactly 
one zero of Fn(x), in agreement with Fig. [1] and Figs. 
1,2 of Ref. [23 ■ Hence, F„(x) has its zeros and poles 
interlaced on the real axis. That would also conclude any 
numerical method of computing F{x) through Eq. dTTI) . 
because of an unavoidable cutoff at some n = N ^ 1. 

The above conclusions remain valid also in the limit 
ri — !> oo. A point of crucial importance is that the in- 
equality (1121) survives the limit as the sharp inequality 



(19) 



for all fc > 1. The sequence in Eq. ([TT|) converges to a 
Mittag-Lejfler PFD, 



F{z)=ao + J2 



fc=i 



6' 



(20) 



defining a meromorphic function in the complex plane C 
with real simple poles and positive residues 



0<Mk 



E 

.;=o 



{1 + 1)1 



(21) 



The series is absolutely and uniformly convergent in any 
finite domain having a finite distance from the simple 
poles , and it defines there a holomorphic function of z 



[z € C here and below has no relation to z in Eq. ^] 
One obtains at once that 



dF{x) 
dx 



< 0. 



(22) 



Because F{x) monotonically decreases from +oo to —oo 
between any its subsequent poles £,k and ^fc+i, there is ex- 
actly one zero of F{x) for x € {£,k, £.k+i)- As a byproduct, 
the spectrum in each parity eigenspace B± is necessar- 
ily nondegenerate. Eventually, the knowledge of another 
OPS, {-Pn~^^(a;)}5^i, enables one to determine the ex- 
pansion coefficients of a physical state described by dH) 



(^n^x) 



Pt^\x) 



(23) 



We prove the above results in the forthcoming section 
im which is divided into two subsections. First, subsec- 
tion lll Al deal with the case of an arbitrary large but finite 
n = N. The limit A'^ — cx) is then considered in subsec- 
tion III Bl Section mil illustrates some of our findings on 
the exactly solvable case of the displaced harmonic os- 
cillator. In Sec. IIVI our results are then extensively dis- 
cussed from various angles. Subsection IIV Al gives a com- 
parison of the properties of our F with those of Braak's 
functions G±. In subsections IIV Bl and IIV Cl recent claims 
re gard ing solvability and integrability of the Rabi model 
[3ril2| are critically examined. Compatibility of our re- 
sults with some other results of the theory of infinite con- 
tinued fractions and complex analysis is demonstrated in 
subsection IIVDI Subsection IIV El gives an overview of 
some open problems. We then conclude with Sec. |Vl 
Some additional technical remarks are relegated to Ap- 
pendix |X1 



II. PROOF OF THE MAIN RESULTS 

According to the Wallis formulas (Eqs. (III. 2.1) of Ref. 
[281 ]: Eqs. (4.2-3) of Ref. [H]), given a three-term recur- 
rence ([3]), the infinite continued fraction in Eq. (jS]) can 
be expressed as the limit 



ro = lim — , 



(24) 



where the nth partial numerator An and the nth partial 
denominator Bn are determined as linearly independent 
solutions of the recurrence 



A„ 



O.nAn-1 
O.nBn-1 



bnAn-2, 
bnBn^2i 



(25) 
(26) 



where n > 1. The A^s and Bn's are differentiated by 
the initial conditions: 



= 1, Aq^ 0, B_i = 0, 



Bn 



1. 



(27) 
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In an intriguing and peculiar world of infinite continued 
fractions, the respective recurrences (PS)) and (|26p are 
essentially identical to the initial three-term recurrence 
([3]) (up to the change a„ — >■ — a„ and the omission of the 
n = Q term). For the Rabi model we have (Eq. (37) of 
Ref. [23], or Eq. ([10]) herein above) 



1 



{n + 1) \k 



bn = 



1 



n + 1' 



c„ = -[n±(-l)"A]. 

K 



(28) 
(29) 



A. Arbitrary large but finite A'^ 

In order to prove the properties of Fn{x) defined by Eq. 
(ITT]) for an arbitrary n, together with the properties listed 
below, it is sufficient to prove that each of the recurrences 
([TU)). (^5)) . (1^ can be transformed into a recurrence of 
the type 

P„(x) = (a; - Cn)Pn-l{x) - XnPn-2ix), (30) 

P_i(x)=0, Poix) = l, (31) 

where n > 1, the coefficients c„ and A„ are real and 
independent of x, and A„ > 0. Obviously, the above re- 
currence defines a family of polynomials {Pn}^=o with 
degree Pn = n. According to Favard-Shohat-Natanson 
theorems (given as Theorems 1-4.1 and 1-4.4 of Ref. [28|'). 
satisfying the above recurrence is a necessary and suffi- 
cient condition that there exists a unique positive definite 
moment functional £, such that for the family of polyno- 
mials {Pn} holds 

£[1] = Ai, £[P^{x)Pn{x)] = A1A2 . . . A„+iJ„„, (32) 

m,n — 0,1,2,... and 6„in is the Kronecker symbol. 
Thereby the polynomials {Pn} form an OPS [30] • Be- 
cause A„ > 0, the norm of the polynomials P„ is positive 
definite, C[P^{x)] > 0, and C is positive definite moment 
functional (p. 16 of Ref. (H). 

With a„ and &„ as in Eqs. (^5)) . a substitution _B„ = 
(— l)"P„/(n -I- 1)! transforms the three-term recurrence 
([^5]) into the recurrence of the type ([5D|) and (I5T]) with 



— , Xn — Xn 



(33) 



where c„ has been defined by Eq. ([29]) and x = e/ k = 
E/g. A similar substitution A„ = (-l)"S'„/(n + 1)! 
transforms (pS]) into the recurrence 

S'„(x) = (a; - c„)5'„_i - A„S'„_2 (34) 

with A„ = n, but with a "wrong" initial condition 

S-i = -1, 5o = 0. (35) 

The latter is not of the required type pil) . Note that a 



recurrence of the type (p4l) yields 
5*1 — Ai, 

5*2 = (x - C2)Ai, 

5*3 = (X " C3)(x - C2)Ai + A3A1, 

S3 = {x - C4)(x - C3)(x - C2)Ai + {x ~ C4)A3Ai 
+ {x - C2)A4Al, 

(36) 

Therefore, a further substitution Sn = XiQn-i trans- 
forms the recurrence (l34l) into 



Qnix) = (x - c„+i)(3„_i(x) - A„+iQ„_2(a;), (37) 

where n > 1, with the "correct" initial conditions 

Q-i=0, Qo = l. (38) 

The recurrence ([57]) together with the initial conditions 
is now of the type ([50)1 and pip with 



Cn+l, 



Xn = A„+i = n -I- 1. 



(39) 



Eventually, a substitution (/>„ — > (/)„/n!, transforms the 
initial recurrence (1101) into 



^^^+i = (^-c„)<^^-A„0±_i. (40) 



The recurrence for (f>^ is again of the type (|30|) and (|3T 
with 



Cn-l, 



A„ = A„_i = n - 1, 



(41) 



where we set A 1 = Aq = 1 7^ for n = 1. Note in 
passing that Ai enters the recurrence (|30|) only in the 
product AiP_i, where P_i satisfies the initial condition 
([5T|) . Therefore we have the freedom to set Ai at our will. 
The respective recurrences for 0„'s, P„'s, and Q„'s have 
all been shown to be of the type 

P(")(x) = (x - c„+„)p(^\(x) - A„+„Pil(x), (42) 
P[i>{x)=0, Pt\x)^l, (43) 

where the coefficients Cn and Xn+a are real and indepen- 
dent of X, and Xn+a > for n > 1. One has a = —1, 0, 1 
for (/)„'s, P„'s, and QnS, respectively. 

We continue to denote the polynomials of the OPS for 
a = by {PnjJ^Lo- They determine the denominators 
P„'s in Eq. The respective monic OPS with a = 

— 1,1 are called associated to Pn's and will be denoted by 
{^n"^}^=o (see Sec. III-4 of Ref. [28|). Because 

A..<Zl£££l. B^.tir^^ (44, 



(n + iy. 



[n + iy. 



it follows at once that the ratio (|24l) can be expressed as 
the limit of the ratios of the orthogonal monic polynomi- 
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als 



ipn{x) has exactly n points of increase, a;„fc, 



To = lim 



Pn{x) 



(45) 



The properties hsted below Eq. pT|) follow straightfor- 
wardly from the classic theory of orthogonal polynomials 
(see esp. Sees. 1.4-6 and III. 1-4 of Ref. Hi]). The ze- 
ros of the polynomials of any OPS are real and simple 
(Theorem T5.2 of Ref. Furthermore, the zeros of 

any two subsequent polynomials Pn(x) and Pn+i{x) of 
an OPS mutually separate each other (Theorem 1-5.3 of 
Ref. 28]). The separation property of zeros ([T^ follows 
from Theorem III-4.1 of Ref. [H. Eqs. pi|and (fTC]) 
follow from Eqs. (1-5.6) and (1-6.2) of Ref~l28(, where 
we have assumed C[V\ = /zq = 1. The partial fraction 
decomposition (ITSI) follows from Theorem III-4.3 of Ref. 
[28j . The positivity of M^k in Eq. dTT]) follows from the 
Christoffel-Darboux identity (Eq. (1-4.13) of [28]), 

P'^+,{x)P,,{x) - P:,{x)Pn+i{x) > 0, (46) 

which for x ~ Xnk reduces to 

Pn+l{Xnk)Pn{Xnk) <0, (47) 

where the prime denotes derivatives. The 2nd of Eqs. 
(IT71) follows from Theorem 1-4.6 of Ref. Thereby 
our results for any finite n — N have been proved. 

Note in passing that the coefficients M^^'s in the PFD 
(1151) satisfy, in virtue of the Gauss quadrature formula 
(Eq. (II-3.1) of Ref. HI), 

n 

J2Mnkxlk=l^i (?-0,l,2,...,2n-l), (48) 
fc=i 

where /i;'s are the corresponding moments of the positive 
definite moment functional, C[x^] — fii. (The positivity 
of C implies fi2i > 0, but not necessarily /i2;+i > 0.) 



B. The limit iV oo 

According to the representation theorem (Theorem II- 
3.1 of Ref. \2§\), the weight function -0 of the positive 
moment functional £ (also called distribution function 



Clx'^] 



(n = 0,l,...), (49) 



is the limit of a sequence of bounded, right continuous, 
nondecreasing step functions ?/;„(a;)'s, 

V'n(a;)=0 (-00 < X < a;„i), 

tpnix) = M„i + . . . + Mnp {Xnp <X< Xn,p+l), 
^nix) ^ Ho {X > Xnn)- (50) 

Consequently 



• the discontinuity of ipn (x) at each Xnk equals Mnk 
(fc = l,2,...,n), 

• at least the first (2n — 1) moments of the weight 
function ipnix) are identical with those of tp{x), i.e., 

dMx) = Mi (? = 0, 1, 2, . . . , 2n - 1). (51) 



Obviously, for any z £ C different from the zeros x„fc's 
the PFD in Eq. (IT5t can be expressed as 



Pn{z) 



E 

k=l 



M„k 



d-tpnix) 



Z — X 



(52) 



According to Hamburger's Theorem XIF [31], the func- 
tion 



dip{x) 



(53) 



where the Stieltjes integral measure dip has been defined 
through the limit of ■!/'„(a;)'s, is a regular analytic func- 
tion in any closed finite region fl of the complex plane 
which does not contain any part of the real axis. The 
convergents in Eq. (j52p converge uniformly to f{z) in 
ri. According to Definition III- 1.1 of [28l], the infinite 
continued fraction in Eqs. ([8|) and (|24|) then converges 
and 



F{z) ^ao + 



dip{x) 



(54) 



So far we have mostly summarized the relevant clas- 
sical results of Hamburger j3l|. A point of crucial im- 
portance in our case is that the resulting Stieltjes mea- 
sure dip{x) = ipix) — ip{x — 0) is necessarily discrete. 
(Here ip{x — 0) denotes the left-side limit of V' at x, 
ip{x — 0) = limx^^x,x„<x'>P{xn)-) To this end, we first 
show that the set of zeros Xnk extends beyond any bound 
up to at +00. Denote 

a = lim ^j, (55) 



where ^j's are the limit zero points defined by Ec^. (|14l) . 
According to Eq. (IV-3.7) of Ref. [11], a sufficient con- 
dition for cr = cx) is that 



lim c„ = cxD and lim sup 



An+l ^ 1 



< 



(56) 



In the present case of the Rabi model, with c„ and A„ 
defined by Eq. ([55]) . one has A„+i/(c„c„+i) — 0{n^^). 
The conditions (1561) are then obviously satisfied. 



lim Cn — 



lim sup 



An+l 



n— s-oo CnC-n+l 



(57) 
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Now the condition cr = oo ensures that the limit zero 
points defined by dH]) are ah distinct, i.e., Eq. 
holds. Indeed, if £_k — £.k+i for some fc, then is a limit 
point of ^;'s (Theorem II-4.4 of Ref. According to 

Theorem n-4.6 of Ref. [28i], if = Cfe+i for some /c > 1, 
then 



^fc = a = lim ^j. 



(58) 



(Such a separation of the limit zero points ^^"■''s applies 
also to the other two OPS with a ^ 0.) Because of 
the sharp inequality (fTO]). the Stieltjes weight function ij; 
satisfies ip =const on any interval x G (^fc,^fe_|_i). Con- 
sequently dip is a discrete measure. Moreover, the mea- 
sure dip is unambiguously determined (see footnote 50 
on p. 268 of Ref. |31[). The determinacy of dip, and the 
Stieltjes weight function ip (assuming the normalization 
ip{—oo) — 0), follows also independently from Carleman's 
criterion which says that the moment problem is deter- 
mined if (cf. Eq. (VI-1.14) of Ref. [H) 



1=1 



-1/2 

= OO. 



(59) 




Figure 2. Analyticity domain of F{z) guaranteed by Hambur- 
ger's theorem [3l[ extends to the half-planes above and below 
the real axis at the distance of at least S from the latter as 
showed by the dashed lines. In the case of the Rabi model 
there is possible to extend the domain of analyticity of F{z) 
through the real axis in any interval {^k, Cfc+i) and to connect 
the upper and lower half-planes. 



The latter is obviously satisfied in our case. 

Now the support of the measure induced by ip, or 
briefly the spectrum of ip, is precisely the set of all ^fe's 
(pp. 63 and 113 of Ref. [2^). Therefore, by the very 
definition of the spectrum of ip (p. 51 of Ref. [1^), the 
residues in ([51]) are strictly positive, < dip{^k) = 
which proves Eq. (I?!]) for z G having a finite distance 
(5 > from the real axis. As the result, the Stieltjes 
integrals in Eqs. ([55]) and ((M)) reduce to infinite sums. 

In what follows we show by the Stieltjes- Vitali theorem 
(cf. p. 121 of Ref. jH; p. 144 of Ref. [13) that the 
convergents in Eq. (|52p converge uniformly to a regular 
analytic function f{z) represented by Eq. (|53p not only 
in any closed finite 57 n M — 0, such as in Hamburger's 
Theorem XIF |31J, but also on the real axis for z = x £ 
(^fc -|- 26,£,k+i — S), where 5 > is sufficiently smah real 
number (see Fig. [2]). To this end we remind that the sum 
over all M„fe's satisfies Eq. p^ . Therefore (cf. Eq. (56) 
of Ref. [13]), 



fc=i 



1 

^ A' 



(60) 



where 6 = min | z — Xnk \ ■ Now select any pair of subse- 
quent limit zero points ^fc and £,k+i- Because the limit 
zero points are separated, there exists some infinitesi- 
mal (5 > so that S^k + 25 < S,k+i — 5. Now consider 
z = X G {£,k + 2S,^k+i — 5). Because each sequence Xnk 
is decreasing with increasing n, Xnk ^ i£,k + ^i^k+i) for 
sufficiently large n > N. Then min |a; — Xnk\ > 6 > 0. 
This establishes a uniform bound on the convergents for 
z — X e (Cfe + 26,£_k+i — <^)- The latter bound can be 
obviously extended to any closed region in the complex 



plane bounded with a semicircle of radius S at + S and 
a semicircle of radius 6 at ^k+i and having a distance at 
least 6 from the real axis for Re z G (— oo, -I- 6] and Re 
z S [^fc+i, oo) (see Fig. Analyticity is then established 
by the Stieltjes- Vitah theorem \2^,^. Thus F{z) takes 
on the form of the Mittag-Leffler PFD which is ab- 
solutely and uniformly convergent in any finite domain 
having a finite distance from the simple poles ^^'s. 



III. EXAMPLE OF THE DISPLACED 
HARMONIC OSCILLATOR 



The recurrence ((72)) has for A = 0, i.e., in the case of 
the displaced harmonic oscillator, the solution [3] 



(61) 



where ( — e + k'^ and are associated Laguerre poly- 
nomials [13 (note different sign of k compared to Eq. 
(2.16) of Schweber [|). Therefore, 

P^-i)(a;) = n!K-"4'^"+'='-")(K2). (62) 

Let instead of the recurrence pop polynomials Qn satisfy 

Qn{x) = (x - — -] Qn-lix) - ^ Qn-2{x). (63) 

\ P J P 



Then 



Qn{x)^p-''Pn{pX + q) (p^O). (64) 



7 



Given Eqs. ((62|) . (I63l). and (l64l). 

Pn{x) - - (I/k)] = n\K~"4r-"H^^)- (65) 

The generalized Laguerre polynomial of degree n is 



J=0 



(66) 



which implies 



J=0 



(67) 



where the (generalized) binomial coefficients 

Eq. (|67p shows that is a sum of polynomials in ^, 

0n = E(-l)"-^'7;7^n(C-O, (69) 

where C = I corresponds to the eigenvalue 

e; = ? - ^2 (70) 

of the displaced harmonic oscillator 2]. Each 0„ is also 
a polynomial of degree n in energy parameter x, which 
can be verified for first few coefficients explicitly: 



h 

t>2 = K 



2 , C(C - 1) 



a; 
2! 















-3!^ 


2 



C(C-i)«' , C(C-i)(C-2) 
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a; 






X 


1 






H) 






3! 






~ 2! " 


3k 



(71) 



Note that (/)i/0o = 2;, which is exactly the n — Q part of 

(E 



IV. DISCUSSION 

Our recent work has been driven by the curiosity as 
to what extent the recurrence coefficients a„ and &„ of 
a model from TZ determine the model basic properties 
[23I [25} . Earlier we looked at the problem from the per- 



100 



^ 



-100 



Figure 3. Plot of ^'(C) corresponding to the recurrence (|72p 
for g = 0.7, A = 0.4, and a; = 1, i.e. the same parameters as 
for G±{C,) in Fig. 1 of Ref. shows corresponding zeros at 
w -0.217805, 6.29563 x 10"^ 0.86095, 1.1636, 1.85076, etc. 



spective of the minimal solutions [20| of the recurrence 
([3]) and showed that the spectrum of any quantum model 
from TZ can be obtained as zeros of a transcendental func- 
tion F [Hjll^li^. In the present work we took a comple- 
mentary view and analyzed in detail the analytic struc- 
ture of Schweber's quantization condition '5]. We showed 
that the function F originally defined by Eqs. ^ and 
([7]) can alternatively be represented by the Mittag-Lejfler 
PFD (PU]) with repelling zeros and with positive residues 
Aik defined by (PT|) . The latter enabled to prove the 
monotonicity of -^(2:) and nondegeneracy of the spectrum 
of the Rabi model in each parity eigenspace B± . The non- 
degeneracy is not knew result - it was initially obtained 
by Kus |19j within the framework of Frobenius's analysis 
of regular singular points [l^l . Yet it is both stimulating 
and inspiring that our approach based on OPS's enabled 
us to prove the basic analytic property of the Rabi model 
independently, without any recourse to its earlier proof. 

We have presented our results while treating the spe- 
cial case of the Rabi model. However, as obvious from the 
proof, our results remain valid for any model of the class 
TZ which can be reduced to the recurrence of the type ([50)1 
and pil) and which satisfies the conditions (more 
general sufficient conditions than (j56p can be found in 
Sec. IV-3 of Ref. [1^.) Essential for obtaining our re- 
sults was to start with the pair of the recently uncovered 
parity resolved three-term recurrences PU)) (cf. Eq. (37) 
of Ref. |23| ) . The recurrences are different from the orig- 
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inal recurrence for the Rabi model 

/n(C) ^ 1 ^ 



(n+1) 
where /„(() is given by 



71+1 



= 0, 



1 / A2 
2k \ n — 



(72) 



(73) 



where k — g/uj and A = ji/uj are as in Eq. pH)) (cf. Eq. 
(A8) of Schweber [2] , which has mistyped sign in front of 
his bn-i, and Eqs. (4) and (5) of [3]). The dimensionless 



is the same as in Sec. 



energy parameter <^ 

mi 

Because /„ in (j73p contains ^ both in a numerator and 
in a denominator, the recurrence (|72l) does not reduce to 
that of the type ^ and (EH) obeyed by OPS's. Never- 
theless, as shown in Fig. [Sj the transcendental function 
F{Q = —/(C) + ro corresponding to ([72]) still displays a 
series of discontinuous branches monotonically extending 
between +00 and —00 and the spectrum can be obtained 
as zeros of F(C) (cf. Fig. 1 of Ref. [Ull^). However, 
one can no longer guarantee that all the residues Mk 
are positive nor ensure the sharp inequality (jl9l) . Note 
that a kind of a Mittag-Lejfler PFD (PO)) is rather gene- 
ral (cf. Eq. (67) of Grommer 132]). It also applies to 
OPS's which distribution function has a bounded denu- 
merable spectrum with a finite number of limit points 
(cf. Theorems 3.2 and 5.4 of Ref. (H). 



A. A comparison with Braak's functions 

In a recent letter [sj , Braak claimed to have solved the 
Rabi model analytically (see also Viewpoint by Solano 
[12]). He suggested Q that a regular spectrum of the 
Rabi model was given by the zeros of transcendental func- 
tions G±{(), 



G±(c) = E^»(c) 



It 



(74) 



Here the coefficients Kn{C) were obtained recursively by 
solving the Poincare difference equation (j72p upwardly 
starting from the initial condition 



= — ao = 2k 



1 

2k 



c 



c 



(75) 



Braak argued that between subsequent poles of the term 
in the square bracket in ([71)) a.t ( = n and C ~ n + 1 the 
function G'±(C) takes on zero value: 

• once - by implicitly presuming that at one of the 
poles G±{Q goes to -|-oo and at the neighboring 
pole goes to —00, with a monotonic behavior from 
+00 to — cxD between the poles; 



• twice - implicitly presuming that G± (C) goes to one 
of ±00 at both subsequent poles of the term in the 
square bracket in (fMl) . and in between the poles it 
has rather featureless behavior, e.g., similar to the 
cord hanging on two posts; 

• none - occurs under the similar circumstances as 
described in the previous item, if the "cord is too 
short", e.g., it does not stretch sufficiently up or 
down as to cross the abscissa. 

However, the above behavior can be merely regarded as 
an unproven hypothesis. There is no proof in Ref. Q 
that the above behavior is the only possible behavior. In 
this regard, Eqs. (|7^ and ([75)) show that Kn{C) is a com- 
plicated function given as a sum of polynomial fractions 
with increasing polynomial degree in both the numera- 
tor and denominator. Therefore, each i^n(C) has its own 
zeros and poles structure. Consequently, it cannot be ex- 
cluded that between any two subsequent poles of the term 
in the square bracket in (TMl) the function G±(C) would 
display much more complicated behavior than that as- 
sumed by Braak [3]. The same objections apply to the 
proposed functional form of G±{C) in Eq. (6) of Ref. [1], 



G±(C)-G"±(C) + 



71=0 ^ 



(76) 



where Gj_(C) is entire in (. There is no proof in Ref. Q 
that Kn{C) 7^ nor that Kn{C) is free of poles [cf. Eqs. 
dZD) and dlSD]. Actually the Mittag-Lefiler PFD (gnj) of 
our F function is much stronger result than the unproven 
Eq. jZe)) of Braak Q. Whereas Braak [3| cannot say 
anything about the residues h^, the residues A^^ in the 
Mittag-Lcfifcr PFD (1^ are all positive and given by Eq. 
()2ip . There is also no entire function contribution in our 
PFD ([^0)) . Furthermore, the F function has exactly one 
zero between any two of its subsequent poles. 



B. Solvability 

A typical quasi-exactly solvable Schrodinger operator 
can be expressed as a polynomial of degree at most two 
in the generators of the sZ(2,M) algebra [l^ [35j, 



(77) 



where qkrm Qrn^ a-i'G some real constants, 



J!!=a„ J^' = z9,--, J^^z^d,-nz, (78) 

n is an integer, and 

[Jl,r]=2j:^, [Jo",J^] = ±J2. (79) 

The Rabi model allows for such a representation in 
terms of the generators of s/(2,R) algebra only for the 
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energies satisfying Eq. ^ (cf. Eq. (12) of Ref. [17|), 
which correspond to the eigenvalues of the displaced har- 



monic oscillator. The lines / 



are exactly the base- 



lines on which the Juddian exact isolated solutions occur 
p^| - |20j . The latter correspond to the degenerate solu- 
tions of the Rabi model [1^ and coincide with the inter- 
section points of different parity solutions as a function 
of energy 19]. The baselines are the only places where 
such a degeneracy is possible [l^. This shows that the 
very existence of the Juddian exact isolated solutions is 
a consequence of that the Rabi model is quasi-exactly 
solvable. The quasi-exactly solvable eigenstates are ex- 
actly the parity degenerate Juddian isolated exact solu- 
tions [iiSi. 

However, not any quasi-exactly solvable Schrodinger 
operator can be expressed in terms of the quadratic el- 
ements of an enveloping sl(2,R) algebra [S^. Therefore 
one might argue that an exact solvability for other energy 
values would still be possible and the claims of Braak Q 
could be justified. Nevertheless, the possibility of polyno- 
mial solutions can be excluded by recent result by Zhang 
[37| . Indeed, let us consider the differential equation 



X{z)— + Y{z)— + Z{z) 



dz 



dz 



$(z) = 0, (80) 



where X{z) = Y.Skz\ Y{z) = j:tkz\ Z{z) = J2vkZ^ 
are polynomials of degree at most 4, 3, 2, respectively. 
Zhang [231 found all polynomials Z{z) such that Eq. (|80|) 
has polynomial solutions $(z) = nj=i(^ ~ ^j) of degree 
I with distinct roots Zj. Theorem 1.1 of Zhang ^3] yields 
an algebraic conditions on each of the coefficients Vk, k — 
0, 1, 2 of Z{z) in terms of the coefficients X{z), Y{z), and 
the roots Zj. In the case of the Schweber's equation for 
the Rabi model [Eq. (3.23) of Ref. 0), and upon taking 
into account that Schweber's n is twice of ours. 



z{z — 2k) 



+ [2« - k) + (1 - 2C + 4k^)z - 2kz^] '^'^^''^ 



dz 



2kC(1 - z)] <^>{z) = 0. (81) 



Because for the Rabi model Zhang's S4 — — — 0, 
Zhang's condition (1.8) on V2 reduces to V2 = 0, and 
Zhang's condition (1.9) on vi reduces to vi = —lt2, or. 



(82) 



The latter leads t o (iTOll. i.e.. again to the Juddian exact 
isolated solutions |18l-l2(]|. One arrives at the same con- 
clusion if Zhang's conditions [STf] are applied to Eq. (5) 
of [13], which is another variant of (ISTI) . Therefore, the 
Rabi model has no other polynomial solution than the 
Juddian exact isolated solutions. 

A possible signature of exact solutions could be seen in 
the case of the displaced harmonic oscillator. There the 
solution of the recurrence for (/>„ was given by Eq. (1691) . 



In general, rapid growths of the product 11^=0 (C ~ Cfe) 
in dnH) essentially cancels out the 1/n! prefactor, leads 
to a finite radius of convergence of the series for $(2:) 
in Eq. and prevents $(2:) from being an element 

of B. The points of the spectra are characterized by a 
sudden collapse of the degree of , which is in general 
the polynomial of degree n in energy, to a polynomial 
of merely the {I — l)th order for any n > I at the Ith 
spectral point (including / = 0). Indeed, ( = e + k'^ is 
a nonnegative integer I at the points of the spectra and 
the sum over j in Eq. (j67p runs only between j — and 
j = I — I for n > I. Otherwise the product on the right- 
hand side of Eq. ((67| vanishes. The leading (Z — l)th 
order in C is then rapidly decreasing with increasing n as 



(-1) 



n+l-l 



Ik' 



{n + l-l)V 



(83) 



which implies $ G ;B 0. Thus the very same product 

terms nfe=o(C ~ Ck), which initially prevented $(2) from 
being an element of B, come later on to the rescue and 
ensure that $(z) G S for the spectral points C = = ^• 

It appears plausible that (pn's can be expressed as a 
sum of such product terms over the spectrum of a model 
and the points of the spectra would again be character- 
ized by such a sudden collapse of the degree of </>„ also 
in general case (e.g., for the Rabi model). However, even 
in the case of the displaced harmonic oscillator, the (/)„'s 
will be generated by the recurrence (j40p as polynomials 
in X. The leading order of 0„ in x is then provided by 
the polynomial term 



-. n — 1 

sR 



X — — 
nl ^ ^ \ K 

1=0 



(84) 



Starting from the recurrence (|40p. it is highly nontrivial 
to arrive at Eq. (p7|) then. The latter is established 
neither here nor by Braak Q. 



C. Integrability 

In his work on the Rabi model Braak Q suggested 
the following criterion of quantum integrability: If 
each eigenstate of a quantum system with fi discrete 
and /2 continuous degrees of freedom (d.o.f.) can be 
uniquely labeled by /i -I- /2 = / quantum numbers 
{di, . . . , df-^ , ci, . . . , cy-j}, such that the dj can take on 
dim('Hj) different values, where Hj is the state space of 
the jth discrete d.o.f. and the Ck range from to infinity, 
then this system is quantum integrable. 

The Rabi model has /i = /2 = 1 and all degenera- 
cies have been known to take place only between levels 
of states with different parity, whereas within the par- 
ity subspaces no level crossings occur [l^, [S^]. Conse- 
quently, the spectral graph of the Rabi model consists 
of two "ladders" of level lines \p, m) . Each ladder cor- 
responds to an invariant subspace of the Z2 symmetry 
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characterized by p = ±, the parity eigenvalue [iol . |20| . 
Because the global label (valid for all values of g) of the 
Rabi model is two dimensional as /i + /2 — 2, Braak 
3 then concluded that the Rabi model belongs to the 
class of integrable systems. However, the above reason- 
ing is not conclusive. First, the very notion of quantum 
integrability is the subject of ongoing dispute [3^. Sec- 
ond, the nearest-neighbor distribution of energy levels is 
customarily used to distinguish between integrable mod- 
els and chaotic systems |39l - t41j . If a normalized spacing 
between two energy levels is s, then, according to the 
Berry- Tabor universality conjecture (39l . l4l| . integrable 
models are generically characterized by a Poisson dis- 
tribution V{s) = e~*, which peaks up at the zero level 
spacing, ^(O) = 1. On the other hand, according to 
the Bohigas-Giannoni-Schmit conjecture ^4ft]jfor chaotic 
systems is symptomatic a level revulsionjA^. l4l| . i.e., 
V{0) = 0. For instance, for time-invariant systems the 
normalized distribution of level spacings of chaotic sys- 
tems is well approximated by the Wigner distribution 
V{s) = (7rs/2) exp{— 7rs^/4} of eigenvalues of the Gaus- 
sian orthogonal ensemble of random matrices [i^ l4l| . 
Given our result that the limit zeros points repel each 
other [cf. Eq. (IT9|) ]. and that the zeros of F{x) are found 
typically midway between subsequent pairs and ^fc+i, 
one finds support rather for the level repulsion, and hence 
chaotic dynamics of the Rabi model in each of the par- 
ity invariant subspaces B±. Obviously, further work is 
needed to resolve the integrability issue. In this regard, 
our approach based on the properties of a characteristic 
function F{x) defined by Eqs. ^ and ([7]) is well posi- 
tioned to deal with that task. The reason is that, with 
the exception of the very first zero, the presence of each 
subsequent zero of F(x) is announced in advance by a 
preceding pole of F{x) (see Fig. [T]and Figs. 1,2 of Ref. 
[23|). This makes the search of all the zeros of F(x), 
which correspond to the spectrum of the Rabi model, 
quite straightforward. 



D. Compatibility 

With increasing n, the PFD (IT5t . ([5^ define a se- 
quence of rational functions with simple real poles and 
positive residues. Kritikos (§ 4 of Ref. [l^) showed that 
if the sequence of such rational functions 



M„ 



k=l 



(85) 



converges uniformly in a proximity of some point zq 6 C, 
then the sequence converges everywhere in the complex 
plane with a possible exception of the real axis. The con- 
vergence is uniform in any bounded region of the com- 
plex plane with a nonzero distance from the real axis. 
The latter is obviously compatible with our main result. 

In general one cannot always guarantee that, such as 
in Eq. (1^ . M„fc Mk = d?/) > in the hmit n — > oo. 



Indeed, Theorem I in §8 of Grommer (32| merely ensures 
that for any m > 

■IpniXnk - 0) < i^n+yniXnk) < ^n{Xnk) 

= ipn{xn,k+i - 0) < (a;„,fe+i), (86) 

where as usual '0n {x — 0) denotes the left-side limit of tpn 
at X. On taking the limit n oo one cannot exclude that 
V'(Cfc) = "0(6+1) for some fc, and hence (iV'(6+i) = 0. 
The above Grommer's theorem appears to be related to 

the fact that, given the sharp inequality ([T9|) . could 
coincide with one of ^fe < S,k+i- Because sjj}'' < Ci+u 
one can only exclude that two subsequent ^^"^ and ^[.^^j^ 
coincide with a single ^ or S,k+i- If Ci^'' = ^fc+ii then 



dip{S,k+i) = in the Mittag-LeffterFFB (|20| . 

On combining the special cases of Eqs. (1-4.12) and 
(III-4.4) of Ref. [281 for x = Xn+i,k and on substituting 
for Pn{xn+i.k) from the former to the latter, one obtains 



Fn (^ri 



-l,k) 



P^+i(a;„+i,fc) 



E 

.1=0 



Pf{Xn+l,k) 
{l + l)\ 



(87) 



On comparing with the right-hand side of Eq. (PT|) one 
finds that at the support of dij: the left-hand side of ((87)) 
has a nonzero limit for n — > cxd. That implies the sharp 
inequality 



meaning that the interlacing property of zeros ()13p for a 
finite n survives the limit n — > oo. Note that the sum 
such as in Eqs. (PT|) and (|87l) enters also the celebrated 
Chebyshev inequalities (cf. Theorem II-5. 5 of Ref. ^28] ) . 



E. Open problems 

The dynamics and long-time evolution of the Rabi 
model is well understood only for rather weak couplings 
K = g/u! < 10~^, where the Rabi model can be reli- 
ably approximated by the JC model 0, @. The lat- 
ter was originally proposed as an exactly-solvable ap- 
proximation to the Rabi model by applying the RWA 
and neglecting rapidly oscillating counterrotating terms 
The JC model provides the basis for, from the Rabi 
model perspective somewhat misleadingly called, strong- 
coupling regime [1, [l^, [U ■ The latter encompasses the 
cavity quantum electrodynamics and associated with it 
vacuum-field Rabi oscillations of atoms, molecules, and 
quantum-dots in a cavity @ . Long-time behavior of vari- 
ous dynamical variables of the JC model can be described 
by analytic approximations [43] and the dynamics shows 
periodic spontaneous collapse and revival of coherence 
[431 13] . With new experiments rapidly approaching the 
limit of the deep strong coupling regime characterized 
by K > 1 [Tl| . the question of major physical relevance 
is that of the dynamics of the Rabi model for k > 0.1 
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11| . A great deal of insight into the dynamics of the 
Rabi model has been gained by Casanova et al. [ll| in 
the limit ujq/uj <^ 1 by means of an expansion in the 
small parameter wq/w. Photon wave nmxiber packets 
were shown to propagate coherently along two indepen- 
dent parity chains of states and, like in the JC model 
[43I Hj] , exhibited a collapse- revival pattern of the sys- 
tem population [ll|. Nevertheless, still only very little is 
known about the dynamics in a very interesting physical 
region of wq /ti-' — 1- The just established link between the 
Rabi model and the OPS's could improve calculations of 
the eigenvalues and eigenfunctions, which are prerequi- 
site for a reliable calculation of the dynamics of the Rabi 
model, and shed light on its long-time evolution for all 
values of the dimensionless coupling k [llj. 

Carleman's criterion states that the Hamburger mo- 
ment problem associated with the positive-definite OPS 
(pO)) and ((3T|) is determined (i.e., the Stieltjes weight func- 
tion is unique), if the condition ([5U|) is satisfied p8| . 
An alternative sufficient condition for the determinacy 
of the classical Hamburger moment problem is that the 
moments (j49p of the positive moment functional C have 
to satisfy (cf. Eq. (3) of Ref. [H; §5 of Ref. [46]) 



nite continued fractions, was shown to be equivalent to a 
meromorphic function F(x) in the complex plane C ex- 
pressed by the partial fraction decomposition (|20|) with 
real simple poles and positive residues Aik defined by 
([21]). One obtains at once that (i) F{x) monotonically 
decreases from +00 to —00 between any two of its sub- 
sequent poles and ^k+i, (ii) there is exactly one zero 
of F{x) for X G (^fc,^fc+i), and (iii) the spectrum corre- 
sponding to the zeros of F{x) is in each parity eigenspace 
B± necessarily nondegenerate. Thereby the calculation 
of spectra was greatly facilitated. Recent claims regard- 
ing solvability and integrability of the Rabi model 3, 12| 
were critically examined (see subsections II V B I and II V C|) . 
Compatibility of our results with some other results of 
the theory of infinite continued fractions and complex 
analysis was explicitly demonstrated. 
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< — n!. 



(89) 



where 9 and A are two real positive constants. Hambur- 
ger's condition (15^ is also a sufficient condition for the 
uniform convergence of the associated infinite continued 
fraction in any closed domain f2 of the complex plane 
which does not contain any part of the real axis (45| . 
The sequence of convergents in Eqs. and (|5^ then 
converges uniformly to the infinite continued fraction in 
Eq. ([8]) , which ensures the convergence of the latter [i^ 
(see also pp. 214-215 of Ref. ^l])- This raises the ques- 
tion regarding a mutual relation between the respective 
asymptotic behaviors of the A„'s and /i,i's. In the special 
case of On = in Eq. (|30|) for all n g N, Bender and 
Milton [43 showed that \n ^ n implies ^„ ~ n! and vice 
versa. It would be interesting to prove rigorously if such 
a relation holds also in more general case of c„ 7^ 0, such 
as for the Rabi model. 

Last but not the least, an interesting question is if the 
ideas presented here could also be extended to systems 
characterized by a three-term difference equation (|3]) with 
periodic coefficients, like in the Hofstadter problem of 
Bloch electrons in rational magnetic fields [43|. 



V. CONCLUSIONS 

On applying the theory of orthogonal polynomials, 
the eigenvalue equation and eigenfunctions of the quasi- 
exactly solvable Rabi model was shown to be determined 
in terms of three systems of monic orthogonal polynomi- 
als. The formal Schweber quantization criterion for an 
energy variable x, originally expressed in terms of infi- 



Appendix A: Technical remarks 

Assuming that the power series X^^^o^"/-^" ^ 
nonzero radius of convergence \z\ > R > 0, i.e., 
much weaker condition than Hamburger's condition (|89p . 
Grommer arrived from an integral representation (cf. Eq. 
(63) of Ref. HI) 



^ d%l>{x) 

jji Z X 



(Al) 



to the Mittag-Leffter PFB (cf. Eq. (66) of Ref. [33) 







(A2) 



Essential to his arguments was that f{z) remained finite 
and different from zero for any nonreal z £ C. However, 
in virtue of 



dx 



-00 (u-xy + 



< 00. 



q 



(A3) 



the very same is also true if the integration range in Eq. 
(|A1[) were, such as in our case, extended to infinity. 

By making use of the Nevanlinna theorem [49] (later 
rediscovered by Sokal [13] in his improvement of Wat- 
son's theorem on Borel summability; for an extension of 
the Nevanlinna-Sokal theorem to differently shaped re- 
gion see Ref. [llj), which was not known to Hamburger 
at the time he wrote his [45| , one can immediately amend 
Hamburger theorem under his item 4 on pp. 33-34 of [45 1 
and prove that Hamburger's function f(z) is asymptoti- 
cally represented by the power series S{z) not merely in 
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an angular domain e < argz < tt 
bitrary small positive infinitesimal, 
tangent to the real axis. 



— e, where e an ar- 
but rather in a disc 



The zeros of P„(a;) can be associated with eigenvalues 
of Jacobi matrices. The latter are derived from the coeffi- 
cients of the recurrence ((30)) upon writing each A„ > as 
a product of two complex conjugate numbers A„ = ZnZ^, 
and assembling them into the Hermitian matrix 





Z2 


. 


. 





■ 


4 




Z3 ■ 


. 
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C3 ■ 


. 














. 


■ 






. 





. 


. 




Cn _ 



(A4) 



Then the eigenvalues of M„ are the zeros of Pn{x). (Ex. 
1-5.7 on p. 30 of Ref. [H]). 



Using the results of Hamburger , one can prove that 

n 
k=l 

converges uniformly toward 

oo 



(A6) 



within any vertical stripe \t\ < R — 5, where i? > is a 
nonzero radius of convergence of the series for W{t) and 
d is some infinitcsimally small number. The continued 
fraction in Eqs. ([5]) and (1241) converges then uniformly 
to the Borel transform 

F{z) = ao+ W{t) e-^'- dt (A7) 

Jo 

in any closed domain of the complex plane with Im z > 
S>0 (cf. § 3 of Ref. [H). 
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